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Abstract 



The aim of this work is to describe the behavior of a device capable to 
generate high frequency (~ THz) acoustic phonons. This device consists 
in a GaAs-AlGaAs double barrier heterostructure that, when an external 
bias is applied, produces a high rate of longitudinal optical LO phonons. 
These LO phonons are confined and they decay by stimulated emission of 
a pair of secondary longitudinal optical (LO) and transversal acoustic (TA) 
phonons. The last ones form an intense beam of coherent acoustic phonons. 
To study this effect, we start from a tight binding Hamiltonian that take into 
account the electron-phonon (e-ph) and phonon-phonon (ph-ph) interactions. 
We calculate the electronic current through the double barrier and we obtain 
a set of five coupled kinetic equations that describes the electron and phonon 
populations. The results obtained here confirm the behavior of the terahertz 
phonon laser, estimated by rougher treatments Q. 



I. INTRODUCTION 

Since the pioneering experiments of Chang, Esaki and Tsu ||, the tunneling processes 
in semiconductor heterostructures have attracted great interest. This has been increasing 
with the development of modern deposition techniques such as metal-organic chemical- vapor 
deposition (MOCVD) and molecular beam epitaxy (MBE) that can provide artificial nanos- 
tructures. Such structures, coupled to electrical leads, can display many unusual behaviors. 

Most of works have focalized on electron transport phenomena through heterostructure 
devices. The double barrier heterostructure (DBH) permits the creation of several kinds of 
electronic devices such as semiconductor laser diodes, ultra-high-frequency oscillators, and 
many others ||. 

One of the most important problems in the study of resonant tunneling in DBH is the e- 
ph interaction. After the works of Goldman, Tsui and Cunningham |4|,[| , it was realized the 
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importance of electron-phonon interaction on the electronic properties of these structures. 
However, little importance has been given to the study of the phonons generated in this 
process, the way they propagate, their decay processes, etc. 

Our purpose is to study the phonon generation in a DBH under the effect of an external 
applied bias. Our device consists in a DBH made of GaAs-AlGaAs. It was designed in such 
a way that the difference Ae between the first excited level e\ and the ground state Eq in the 
well is of the order of the LO phonon energy fwj\, for a small applied bias V. Here, u\ is the 
T-point LO phonon frequency. The Fermi level Ep in the emitter is such that the excited 
level is above it. With a further increase in V, the levels in the well are lowered relative to 
e f . When the ground state comes under the bottom of the conduction band, the current 
is almost suppressed until the level E\ reaches Ep. To continue increasing V, the current 
begins to flow through the excited level, but since Ae remains less than hui, the phonon 
emission is inhibited. For a given bias V the resonant condition Ae rs fiu)i is achieved and 
the electrons begin to decay to the ground state by emitting primary LO phonons. The 
potential profile and the level positions at this resonant condition are shown in figure |TJ 

For an Al concentration greater than 0.25 |J or 0.3 J7|, these phonons are confined inside 
the well (they can be also absorbed by exciting electrons from eo to E\). 

The process described above, acts in parallel with the decay of primary LO phonons 
due to anharmonicity. One product of this decay is a secondary longitudinal optical phonon 
{LO), the other is a transversal acoustic (TA) phonon |8|||. The LO — TA pair is produced 
by stimulated emission. Therefore these TA phonons, for a bias greater than the operation 
threshold, are coherent and form a beam that we call saser by analogy with a laser. The 
coherence of the device presented here was studied previously |10| by developing a formalism 



similar to that employed usually to study the coherence in lasers [0,0 which take into 
account the competition among the several emitted modes. 



In the early 60's [13], generation of coherent phonons was discussed theoretically and 



observed experimentally ]H , |T5 |. Several experiments and applications were suggested flIS 
These phonons can be generated by intense pulses of laser [^7|-pl|. More details were de- 



scribed in the review of Kurz |22| . In a recent review, Merlin p3[ discusses the generation 



of pulses of coherent THz optical phonons by intense femtosecond-pulsed lasers. 

Coherent acoustic phonons can be also obtained with a wavelength given by the period 
of the superlattice where phonons are optically generated (i.e., 500 A) |[24|| . In our case, 
the system generates coherent sound by electronic means with much shorter wavelength. 
The generation of high-frequency monochromatic acoustic waves by laser-induced thermo- 



modulation was reported by Damen [pq] . The phonon frequency ranged from 2 to 4 GHz 
comparing with 2 THz for our device. 



Generation of phonons by stimulated emission was theoretically considered |2(J and ob 



served experimentally. Bron and Grill p7| reported the direct observation of the stimulated 
emission of 24.7 cmT 1 phonons in a three level system of V A+ ions in AI2O3. Prieur et al. 
p8|,|29[| studied the phonon emission in a two level system of a glass. These phonons had a 



wide frequency distribution around uj = 0.34 GHz. Zavtrak and collaborators p0|-|35[| pro- 



posed a device consisting in a dielectric liquid with small particles or gas bubbles working 
in such a way similar to a free-electron laser. That device could produce coherent acous- 
tic phonons by stimulated emission at a frequency uj pa 2 KHz. Stimulated emission of 
phonons in an acoustic cavity was obtained experimentally by Fokker et al. [|55| , |57| . They 
use the metastable Zeeman-split doublet in ruby. These phonons have frequencies of about 
60 MHz. 

At present, several kinds of amplified coherent beams of bosons are produced by stim- 
ulated emission. Besides the traditional visible or near infrared photons (lasers) and mi- 
crowave photons (masers), other coherent photon beams are produced. High energy lasing 
processes like X ray lasers (called also xasers) p8| , |39l and gamma ray lasers (GRASERS) 
PDfl were studied. 

A very exciting field has been recently opened with the development of atoms lasers 
pT|^4]. Several applications of these novel devices, like high precision gravimetry, are 
proposed. Strongly related, it is the study of vibrational amplification by stimulated emission 
of radiation (vaser) 



Kleppner [[|6] claims against that complicated jargon of acronyms and proposes to sim- 
plify to atom laser, X ray laser, etc. 

In summary, a great effort has been made to produce terahertz phonons using other 
mechanisms different of stimulated emission. On the other side, acoustic phonons have been 
produced by stimulated emission but at frequencies lower than 2 THz. The DBH resonator 
proposed here could produce a continuous beam of terahertz coherent phonons. The ultra 
short wavelength of these phonons permits potential applications. 

Acoustic imaging is a very well known procedure, mainly in medicine, where it is exploited 
the fact that different water contents on several tissues produce different attenuations for 
a signal of about 4 MHz. Today are available commercial equipments that can produce 
three-dimensional images of very high resolution. 

Equipments for acoustic microscopy are less known but they are also commercially avail- 
able. The most used technique is Scanning Acoustic Microscopy (SAM) [^7] that can attain 
resolutions higher than 0.4 fim, taken with phonons of 2 GHz. 

Besides the scanning acoustic microscopy, several techniques have been presented in a 
recent review [|48|| : photoacoustic microscopy (PAM), scanning electron acoustic microscopy 
(SEAM) and scanning laser acoustic microscopy (SLAM). The theory of acoustic imaging 
was also reviewed, schematic diagrams of the different microscopes are presented and several 



industrial applications are discussed in reference |48l . 

The equipments mentioned above work at up to 10 GHz, three orders of magnitude 
lower than the frequency of the phonon beam considered in the present paper. Therefore, 
they have a limited theoretical resolution given by its wavelength that is 1000 times greater 
than the generated by a DBH device. 

To make the acoustic nanoscopy possible, it is necessary to have a coherent beam of ultra 
short wavelength phonons. To obtain images of a small structure, one would need waves 
shorter than the size of the details of the system. 

H. Maris |EJ describes the process of acoustic nanoscopy by mean of echoes of ultra-short 
phonon pulses optically generated. They are detected by measuring the variation of optical 



reflectivity due to the scattered phonon pulses. They work with coherent phonons of about 
2 THz like those produced by our device, but they were not produced by a phonon laser. 

Evidently, a laser is not necessary to produce optical images. Nevertheless, lasers are a 
very useful tool for imaging and up today they are indispensable to make holograms. On 
the other hand, many optical experiments that can be made with incoherent light are done 
more easily and more precisely with the aid of a laser. Therefore many experiments that are 
made today with incoherent phonons could be done better with the aid of a phonon laser. 

The phonon laser described in this paper has the wavelength smaller than 25 A (what 



is about 20 times less than for the device described in |24j]). Therefore, it is possible to use 



it to see structures of about 50 nm. A coherent source of X rays of about 20 nm begins 



to be available now |[50|j51[| . Even if it will be possible to make images of structures at the 
submicron scale by using X rays of a wavelength of about 25 A, they would have energies 
of about 0.5 KeV, what could affect the system to be studied. 

The device presented here could be used to get a hologram at nano-scale. To do this, 
it is possible to build in a GaAs an acoustic analogous of a Michelson interferometer. It 
consists of a thin barrier of AlGaAs designed as a semimirror, placed at 45° to split the 
beam and to force it to pass one half through the part of the sample we are interested in and 
the other across a very pure reference GaAs. After that, we can use mirrors (made of thick 
AlGaAs barriers) to superpose them and record the interference pattern with a detector. 
This pattern could be processed in a computer to get a hologram of the sample showing the 
three-dimensional shape of the nanostructure. In a previous paper |]J we described, as an 
alternative, the acoustic version of a Sagnac interferometer. 

Because of the high sensitivity to all types of lattice defects, the phonon beam could 



be used to give information about the inner lattice structure |52 |, |49j and form part of a 
phonon imaging setup as a source |53| . 

We will show that the amplitude of the phonon beam is with a great precision directly 
proportional to the incoming current. Therefore, the saser amplitude could be modulated 
by modulating the amplitude of the current. This will permit the use of the phonon laser as 



an information carrier between circuit components at very short distances. These "phono- 
electronic" devices could work at smaller distances and lower energies than optoelectronic 
devices. 

If we achieve to make a "phonoresist" , i.e., a material with a soft mode close to the 
frequency of our phonon laser, the saser could be used for nanolithography. Due to the 
short wavelength and low energy of the phonon beam, "phonolithography" would be a 
better option than the methods currently in used. 

The rest of this paper is organized as follow. The Hamiltonian which describes our 
system is presented in Sec. II. The solutions of the Hamiltonian are obtained is Sec. III. In 
the Sec. IV, we show the procedure used to obtain the kinetic equations. We devote Sec. V 
to present our results and conclusions. 

II. THE HAMILTONIAN 

To describe our system, we considered a Hamiltonian composed by a single particle part 
and one representing the interactions. The single particle part can be written as 

H Q = H e + Hi + H 2 + H 3 , (1) 

where H e describes the electrons, and Hi, H 2 , H 3 the LO, LO and TA phonons respec- 
tively. 

Assuming the grown direction as the z axis, the Hamiltonians Hi, H2, H 3 for the phonons 
can be written 

Hi = Y,^Ajw, (2) 
qi 

^2 = (^2 - ihfy) 4 qa 6 2q2 , (3) 
q2 

H 3 = J2(hw 3 -ihK 3 )bl s b 3ci3 (4) 

q3 
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where hui\, hu>2, hu>3 are the energies for LO, LO and TA phonons respectively and 
6j q (6j q . ) are the creation (annihilation) operators for phonons with momentum perpen- 
dicular to the z direction. Two imaginary terms (ihn2, ifiK 3 ) were introduced to take into 
account the decay by anharmonicity of the LO phonons and the escape of the TA phonons, 
respectively. 

To describe the electrons, we use a tight binding approach with hoppings v between 
nearest neighbors. Writing the Hamiltonian in a Wannier basis, in terms of the electronic 
creation (annihilation) operators cj(c r ) at site r, we get 

H e = J2 £ r C t C r + V (4<V + 4' C r) (5) 

r <r,r'> 

where r = (I, m, j) with I, m, j = — oo, . . . , oo. As the system has not magnetic behavior 
the subscript a (for spin) is left implicit. 

Our system has translational symmetry in the perpendicular directions (xy). Therefore 
the Hamiltonian can be uncoupled by expanding the operators c r in plane waves in these 
directions 

c r = £c jk e* k - x - (6) 

k 

We can treat the system as a sum over ID Hamiltonians for each wave vector k perpen- 
dicular to the current (z) direction 

H e = 2 { 2 ej-kc] k c jk + v ( c 5k9'k + 4'kCj-k) } (7) 

k j <jj'> 

where = — Av. The energies are measured from the bottom of the conduction 
band of the emitter and = Bj + e^. The energies Ej are chosen to describe the energy 
profile of the DBH. For the sake of simplicity, we will leave implicit the dependence on k. 

For the z direction, we separate the space in three regions: the dispersion region and two 
semi-infinite one-dimensional chains. On the left we will have planes with energy Sj = for 
j < and to the right Ej = V for j > L + 1, L is the length of the DBH. The corresponding 
eigenstates of these two regions are planes waves. If we disconnect the DBH from the left and 
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right chains we get for the dispersion region the profile of an infinite (even not rectangular) 
well as it is shown in figure |2|. 

For the dispersion region, we diagonalize a three-diagonal matrix of order L corresponding 
to the profile showed above getting the eigenvalues e m and the eigenvectors \m) of equation 



H' ez \m) = £ m \m). 



(8) 



Here, 7i' ez is the part of 7i ez that goes from the beginning of the left barrier to the end 
of the right one. Written in the basis of planes this is 



Ei V 

v e 2 v 
v e 3 



\ 



(9) 



v e L 

where ij = Sj — 2v. 

After the diagonalization we get L discrete levels. We label these levels through the 
index m as m = 0, 1', 2' . . .. For the right chain we rename the planes, from L + 1, L + 2, 
. . . to 1, 2, . . .; and for the left one, from 0, —1, —2,. . ., — oo to 1, 2, . . ., oo. Therefore we 
get a new picture where we have two semi-infinite chains, the first one corresponding to the 
planes j = oo, . . ., 1 and the other one corresponding to j = 1, . . ., oo. Among the L levels 
obtained from the diagonalization of the matrix (|9|), only the two states with lowest energies 
will participate significantly in the electronic transport. The others are far over the Fermi 
level and their effects will be considered later. 

To connect the DBH with the left and right chains we compute the matrix elements 



Vjm = (j\Hez\m). 



(10) 



The electronic Hamiltonian is 
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H e = J2 E A C 3 + H £ mC ] m C m + V (c] Cj+1 + c] +1 Cj) + 

(4ci + 4c m ) }. (ii) 

m 

A diagram representing the electronic Hamiltonian (|ll]) is shown in figure |3|. 
As it is well known ||26|| , the dominant dispersion process in polar semiconductors is due 
to the coupling between electrons and LO phonons. 
The interaction Hamiltonian 7i in4 has three terms 

'Hint = 7~Le—e ~l" 7~(-e—ph ~l~ 7~~tph—ph (12) 

where H e -e, 7~£e- P h an d H p h- P h describe the electron-electron, the electron-phonon, and 
the phonon-phonon interactions respectively. 

7i e -e is taken in the Hartree approximation. That means that the energies e m of the 
levels in the well depend on the accumulated charge. 

The electron-phonon interaction Hamiltonian has the form, 

tte-ph = E # kqi ( C O k _ qi Cl'k&I qi + 4'k C k - q > qi )- (13) 
k,qi 

In order to recover the translational symmetry in the xy plane we take the average on 

qi i 



Ke~ph =J29± (4 k Cl'k6l k + cl'kCOk^lk) ( 14 ) 



where g k measure the effective strength of the interaction and was calculated as in [pi 
The phonon-phonon interaction is given by [0] 



qi 



This Hamiltonian describes the decay LO — > LO + TA and its inverse process (recombi- 
nation). 

As in the GaAs the TA phonons are emitted fundamentally in the [111] direction ||, we 
considered here a device that was grown in this direction. The inside barrier walls act as 
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mirrors of a Fabry-Perot interferometer selecting discrete cavity modes. This mirrors also 
serve for a selection of the saser modes: those running in the grown direction are reflected 
several times between the mirrors and stay longer in the well, while all others are lost. 

With the increasing of the applied bias V, from a certain value, called saser threshold, the 
TA phonons associated with the cavity resonant mode, i.e., q% = 0, stimulate strongly the LO 
phonon decay and their number grow sharply. Once this mode is sufficiently populated, the 
probability of stimulated emission for that mode exceeds the sum of the emission probabilities 
for all other modes. Thus, the TA phonons emitted in this mode will slave the other ones, 
giving rise to a coherent acoustic phonon beam. This process is similar to what happens in 
a laser. Due to the momentum conservation qi = q 2 + q3 and because of q 3 = |l(J we get 

q2 = qi 

The parameter 7 qi appearing in (plj) represents the intensity of the interaction and was 
estimated from ||, |55| and [|56|j . 



Finally, the total Hamiltonian is the sum of both terms: the single particle and the 
interaction parts 

H = J2n o (k)+H mt (k) (16) 

k 

where we had made the approximation that the momentum qx (appearing in expressions 
(|), (D) and (|l"5D) varies in the same range as k. We can see at expression (|1^) that the 
Hamiltonian became, in this approximation, fully uncoupled for each k. This simplifies 
strongly the obtention of its solutions. 

III. THE RESOLUTION OF THE HAMILTONIAN 

In order to solve de Hamiltonian described in the former section we introduce the oper- 
ators 

Z/ii yn.2 ynz 

^1,^2,713 _ „ "lk u 2k u 3 / 17 \ 

u j* - c ik /^-r rzn /zrr ^ Li > 



The operators Q" r }^ n2 ^ creates the eigenstates of H (without the hopping terms) 
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\jkn 1 n 2 n 3 ) = O^ na ^\0) (18) 

formed by an electron at plane j and n\ LO phonons, ri2 LO phonons with parallel 
momentum k; and n 3 T A phonons. 

As it was done before, we will leave implicit the dependence on k. 

Using the operators ([17]) and the Hamiltonian ( fT6|) we can calculate their equations of 
motion. In order to simplify the notation we call n for the set of {ni, n^, n^} phonons 

dO? 



ih- 



dt 

By expanding the eigenstates | $) of 7i in the eigenstates of Ti 



o»,n . (19) 



l*> = £Pitfn>, (20) 
and due to the orthogonality of |jn), the amplitudes can be calculated as 

fi={jn\*)=(0\O?\9). (21) 



Making use of the equations turning out from (|19|) and the previous definition for the 
amplitudes, we get for j ^ 1, 1, m (m = 0, 1', 2', 3' . . . L) 

dv 11 

lh M = £ W + Tj + VP > +1 + VP ^ (22) 

and for j = 1,1, m 
dp^ 

lh ~dt = £ ^ + Tj + VlvP ^' + + vp % + ^ ^mPm, (23a) 

m>2' 



i% db = £ ° P ° + 9y/ ™~ lP "' + T ° + ^ + VmP ^ (23b) 

ih lHr = ^ + aV^h+Tpf + TV + v- iv fi + wiP?, (23c) 

ih^§- = e?p? + Tx + wvft, + v 01 p% + vp%+J2 vimPn, (23d) 

dt m>2' 

ih^ = elp»+v lrn p't + v Im p? + r m , m>2' (23e) 



where 
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e™ = Ej + n\hu)\ + n 2 (huj 2 — 1T1K2) + n 3 (hu 3 — i^/t 3 ) , 
n = {ni - l,n 2 ,n 3 } , 
n + = {m + l,n 2 ,n 3 } , 

t 3 - = 7 {VJFjjf + VW+pf} , 

n' = {m - 1, n 2 + 1, n 3 + 1} , 

n" = {ni + 1, n 2 - 1, n 3 - 1} , 
AT- = ni (n 2 + 1) (n 3 + 1) , 
iV + = (ni + 1) n 2 n 3 . 

The amplitudes are time dependent. We transform them 

1 f°° 

p»(t) = — du>t%{u)e**. (24) 
V 2n J-00 

By doing that, we get the following set of equations 

huja" = e"tf + Tj +va" +1 + va"_ 1 , for j ^1,1, m (25a) 
?iufl" = Ejdj + Tj + fa^r + i^ x ,a" + ^Toa" + X! ^m ^ (25b) 

m>2' 

ftc^a" = £o a o + dV™i a i' + To + "yro a T + ^oi a" , (25c) 

hum™, = e"a" + gV^i + l«o + + Ti' + vj v aj + fi'ia", (25d) 

ftwa" = e"a" + I\ + va 2 + uvio" + voi^" + X! ^1™°-") (25e) 

m>2' 

= e= t£ + T m + v Jm a? + w lm a° m > 2'. (25f) 
Here 1^- is the Fourier transform of Tj. 

In order to obtain the solutions of (p5|), we considered only three channels through which 
the electrons can tunnel. Such channels are the only that are between the bottom of the 
conduction band (e = 0) and the Fermi energy e\ at the emitter, near the resonant condition. 
Therefore, they are the only ones that contribute significantly to the current. These channels 
are: one electron through the excited level and no phonons (channel 0, i.e., n = {0,0,0}); 
one electron through the ground level and one LO phonon (channel 100) and one electron 
through the ground level and a pair LO — TA phonons respectively (channel 011). 
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To solve the system above, we drop the imaginary part of ej. This cope with the problems 
caused by the non-Hermitian Hamiltonian. 

We have to note that the amplitudes a° (for j > 1 or j < 1) remain uncoupled (this 
tunneling channel does not mix with any other). On the other hand, the other two channels 
are coupled due to the ph-ph interaction. To uncouple them we propose the following 
transformation for the amplitudes 

af° = -^(at + a 7 ), (26) 

S 011 = ^ (4 - «7) • (27) 
Making the above transformations for the amplitudes a* 00 and a° n we get two new 
uncoupled channels (called channels + and — ). 

Considering only the channels and ±, the system (EH) can be rewritten as 



va]_ x + (ej - huo) oj + va° j+1 = (28a 

vaf_ x + (ef ± 7 - fiw) + va^ +1 = for j 7^ T, m, 1; (28b 
+ ( £ o _ ^ a o + ^ i;a o + ^ = (2gc 



uaj + (4 00 ± 7 - fcw) a± + v- lQ 4 + X! ^T m «m = (28d 

m^O 

Ul'ia?y + Yl "™l fl l + ( £ 1 - huJ ) a l + VCL 2 = ( 28e 

"oi4 + E w -i a m + ( £ i°° ± 7 - H af + vaf = (28f 

m^l' 

vj m a^ + (4°° ± 7 - hu) a± + v ml af = for m ^ 0; (28g 

u Tm a?+(e^-nw)o^ + v m ia? = formal'; (28h 

vj a± + (e J 00 ± 7 - 4 + -J=a$, + v ia± = (28i 

v Ti' a T + ( £ i' ~~ ^) a ?' + ^| a o + ^/f a ° Ul ' ia i = (28j 
The amplitudes and a° m appearing in equations ( p8|) take into account the electronic 
energy levels different from the fundamental and excited one. To eliminate these amplitudes 

(the fundamental and excited levels are the only ones that participate directly in the LO 
phonon generation) we use the equations (|28c|) - (|28b|) . The resultant system is 
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ua°_i + (e°j - hu) a° + va° j+1 = 

vaf_j_ + (.s] 00 ± 7 — hoSj af + vaf +1 = for j ^ T, m, 1; 
ua| + (ef - ^w) 4 + ^Ti' ?' + ^Ti a i = 
va^ + (e^ ± 7 - fiu;) + + ^Ti^ = 

ujjaf + ui'ia?, + (e? - huj) a\ + va° 2 = 

+ v 01 ciq + (if ± 7 — ftu;) af + waf = 

v Ti' a T + ( £ ?' ~~ ^ luJ ) °i' + "^ a ° y^f a ° Vl ' ia i = ^ 
^To^T + ( £ o°° ± 7 - ^) ^0 + "^| a i' + Uo1 a i = 

where 



. . \ " V Tm V ml 



E 



and 



V Tl= E 



" 4 00 =F 7) 

are renormalized hoppings joining sites j = 1,1, and 





= 4+4 








= £ ? + s? 




= £ r + sf 



are renormalized energies at the sites j = 1, 1 for each channel, where 



4= E 



lm 



4=E 



v 2 



m^O 



(»o; - £i? =F 7) 
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A diagram representing equations (^9|) is shown in figure |j. 



To solve the system (|29|) we proceed as follows 



For the channel and ±, in the region j > 1 the solutions of ( |29a| ) and (|29b| ) are 
respectively 



a° = A%e lk ^ (34) 



a ± = A±e lk±z > (35) 

where ^ is the coordinate of the plane j (zj = ja; a is the lattice parameter), Aj>, A T are 
the amplitudes of the transmitted waves, and k and k are the z component of the waves 
vectors that fulfills the dispersion relation 

hu = a® + 2v cos (k a) , (36) 
hu = el 00 + 2v cos (A^a) ± 7 (37) 

We did not consider the terms e~ lfcnZj because we assume that there are no electrons 
coming in from the right of the DBH. 

In the region j < —1 we have two kinds of solutions. The first one is obtained for the 
excited level that is above the bottom of the conduction band, in this case the solution is 

a° = A R e- ikoz i + A^ z i (38) 

where A R and Aj are the amplitudes of the reflected and incident waves and k satisfies 

hu = Sj + 2v cos (k a) . (39) 

The second kind of solution is obtained for the ground state that is below the bottom of 
the conduction band then, in this case, the solutions are evanescent modes with amplitudes 
A\ and A~^ 



a 



J 



A%e K z i (40) 



where the k ± are positive and fulfill the relations 
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huj = Sj 00 + 2v cosh (k+o) ± 



7- 



(41) 



After replacing (0), (0), (0), and (f4Tf ) in (|29|) we get a system of nine equations with 
the following nine unknown quantities 

a i'i a o\ a o i A-Ri A^, A E , A~^, A T and Aj,. 

The solution of this system permit to calculate the electronic current taking into account 
the e-ph and ph-ph interactions and to obtain a set of five kinetic equations as it is shown 
in the next section. 



IV. THE KINETIC EQUATIONS 

To get the kinetic equations that describe the population of electrons and phonons we 
need to obtain the equations of motion for the electron number operators cjcj and the 
phonon number operators b\bi 



d 



ih—(c)c^j= c]cj,H , ih—(b\bi^= b\bi,H 



dt 



(42) 



We define the average populations for electrons in the excited and ground states as 



nj = ^2($\c]cj\$) for j = 0, 1'. 

occ 

From the definition (p0|) of |$), the electrons average populations can be calculated as 



n 



■J Z^i \rj 
n,occ 



(43) 



For the LO, LO, and TA phonons, the average populations are 



17; 



£W&I$>= £ k n i for z = 1,2,3 



(44) 



n,occ,j 



Using the equations (JE3) we get the set of quantum kinetic equations 
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d ^ = G [,_ G o_ qE (45a) 
;/7} ° (:[., c ( > ■ Gf, (45b) 



dt 

= Gf, - Gt (45c) 



where 



dt 

d -^ = G ?o-Dro («d) 

^§ = G E LO - Ex A) (45e) 



^ = |E^Im{K*Pr} 5 (46) 

n,occ 



is the electronic input rate from the emitter to the excited level, 



G °v = \ E^'I^K'P?}- (47) 



^' n i i 



is the output rate from there to the collector. 

Gy = \ E 9 Im {po + >? } v^iTT, (48) 

n,occ 

represents the net balance between emission and absorption of LO phonons via electron 
transitions, 

^ = !£^oIm{po n Vr}, (49) 

n,occ 

is the input rate from the emitter directly to the ground state and 

G 'o = |E^iIm{K>o}- (50) 

n,occ 

stands for the output rate from the ground state to the collector. 

The balance between the decay of one phonon LO by emission of a pair of LO and TA 
phonons, and its inverse process (i.e., recombination) is given by 

Gfo = |E7lm{pf>f }VN, (51) 

'* n,occ 

where N = (rii + 1) (n 2 + 1) (n 3 + 1). 



The decay rate of the LO population turns out to be 

D ro = 2k 2 E {pfPj} n 2 = 2^2^- (52) 
Finally the escape rate of the TA phonons is given by 

E TA = 2k 3 {pTp?} ™3 = 2K 3 nj. (53) 

A diagram for the input and output electron rates is shown in figure [5|. 

The kinetic equations (|4"5]) can be approximated by a system involving only populations. 
To do that we have to replace the amplitudes from ( p2|) and ( p3|) (See the appendix for 
a detailed calculation). Then, we get a new set of kinetic equations 



dt 

~dT 

dn{ 



■- G - R v n^ - R w (54a) 
-R n^ + Rw (54b) 
R w - i? 7 (54c) 



dt 

d ^ = - K2 n- 2 + R i (54d) 
dn^ 

— = -n 3 n 3 + R 1 (54e) 

where 

Rw — w \^h 7 (jh + 1) _ % W] (55) 
Rj = 7o[^i("2 + l) + - {rh + 1)^2 rh)- (56) 

G is the input current, Ry and i?o are the escape rates for one electron through the right 
barrier and w is the transition rate due to LO phonon emission estimated in the appendix 
(and more roughly in [|1]). The other parameters are: the decay rate of a LO phonon j a , 
the decay rate of a LO phonon k 2 and the TA phonon escape rate k 3 . 

It can be seen in the expressions (|5^) and ([56]) that the Bose factors were obtained. That 
means that the phonons are produced in a stimulated regime. For the case of TA phonons, 
the Bose factor % + 1 is thousands times greater than the others. 
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V. RESULTS AND CONCLUSIONS 



The results presented here correspond only to the stationary case, i.e., (t) = aj (u) e lult . 
In this case, all the populations fl43|) and (f44|), and the rates (^)-(^) are time independent 
and they can be written as the same expressions of a" (u) instead of {t). The amplitudes 
a" (a;) are the solutions of the system (|29|). 

The following parameters were used in our calculations. The barriers are 42.4 A wide 
and 300 meV high (38 % of Al concentration). The well has a width of 203.4 A. The Fermi 
energy at the emitter was fixed at 15 meV . The values estimated for the decay rate K2 and 
the escape rate were 30 ps~ l and 0.05 ps~ l respectively. The other parameters used are: 
w, that depends on the applied bias but it is of the order of 2.5 ps^ 1 , 7 Q = 0.11 ps~ x taken 
from ||, and g, that also depends on the applied bias and it is of the order of 3.2 ps~ l . The 
value of 7 was taken equal to 0.65 meV and was calculated from the experimental value of 
7o using f55[ and |56fl . 



The transmittance through the DBH is calculated from 



T 





2 + 




2 

+ 


2 




to 



(57) 



A plot of the transmittance as a function of energy is shown in figure |6|. It can be seen 
that when the ph-ph interaction is considered (taking 7 ^ 0), a central peak appears. The 
other two peaks correspond to the polaronic branches already shown in previous works |57 



Once the transmission probabilities through each channel are obtained, the output elec- 
tronic current (equal to the input current at the stationary regime) through the structure 
can be calculated by summing over the occupied states at the emitter. This sum can be 
easily transformed into the integral 



S m* J (4 - 3 (e) P (e) 



G (V) = -^^r / ^ de (58) 

v ' 2tt 2 h 3 J 1 + R + T K ' 



where S = 0.5 • 10 -3 mm 2 is the device area, m* = 0.067 m e is the effective electron mass, 
R is the reflectance and 
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3{e) 
p{E) 



Av 1 

1 



.4 



sin (k' aj + sin (k + aj + A T sin (k aj 



2v sin (koa) 

Instead of using the expression fl58|) to calculate the electronic current G, it is calculated 
from the sum of the terms (f46|) and ([49|). 

The difference between the currents calculated for 7 ^ and 7 = is appreciable only 
in the dashed region of figure [7|, that it is amplified in the inset. 

Finally, the saser intensity and the TA phonon population are shown in figure [8]. This 
curve is obtained by solving the kinetic equations Q54|). 

In conclusion it was obtained here a set of kinetic equations using a full quantum treat- 
ment. This set is different from those presented in a previous work [lj where the kinetic 
equations were derived phenomenologically. Besides, we consider for the first time all the 
relevant phonons and its interactions in the calculation of the current. 

It was shown that when the phonon-phonon interaction is considered, only small modi- 
fications are observed for the current. 

The accurate full quantum calculations for a DBH phonon laser confirm the results 
obtained using rougher phenomenological methods 
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APPENDIX: 



In this appendix, we will show the way to obtain each term of equations ([54]) from (plop. 
To illustrate the procedure, we will limit us to show only two examples. 
We begin with the term 

g ? = |e w m K^}- (ai) 



n,occ 
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The channel 000 is the only one that contribute significantly to the sum appearing in 
1|) (see Fig. [| and Fig. Taking this into account, the expression ( |A1| ) can be rewriten 



as 



(A2) 



As the amplitudes appearing in ( |A2| ) are time dependent, we substitute one of them by 
their Fourier transform using (|24| ) 



Im {pfpi,} = Im j-^= 



(A3) 



2n j—oq 

To obtain an expression for the amplitude af (a;), we use ( |29ej ) and neglect the term 
containing the renormalized hopping v~ x in it. Taken into account the dispersion relation 
(BSD, (P2cj) and using to calculate the amplitude in the plane j = 2, we get 



af (uj) = E vi af, (uj) 



(A4) 



where 



-l'l = v vi 



E? - ve~ ik ° a 



Sf - 2uE?cos (k' a) +v 



Replacing (|A§ in (|A3|) 

Im = -L Im | y°° H ri a?; (uj) due^p ^. 



(A5) 



The term Eyi inside the integral varies very slowly in the range where the amplitude 
a\, (uj) is not null. Therefore, we can take it out of the integral. Now, performing the integral 
over uj and taking the imaginary part we get 



Sf - 2vS? cos (k' a) + 



sin ( A; n a 



Pi' 



(A6) 



Defining ni' = Y^occ \Pi>\ 2 we get 
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where 



E? 2 - 2vE? cos (k' aj + 



sin (kQa^ . 



When the renormalization efects are not taken into account, i.e., E^ = 0, the expression 



obtained here for Ry is the same as that appearing in [59|. 

The other term to be analyzed here is the rate of LO phonon emission (fJH 



(A7) 



First, we expand the term inside the sum and change the range over which the ri\ runs 
to obtain (|55|). That is 



E Im {pfpy} = £ KPo 1+1)n2H X' V^+T - rfT 1 " 1 ^^) • (A8) 



(ni— l)ri2rig 



L i."i+l)n 2 n 3 j n ^ e |- erm Q f 



Writing p 
ing the san 
the terms that contain the factor g) we get 



in terms of its Fourier transform a, 



(ni+l)ri2nj 
: 



doing the same with p^ 1 1 ^ ri2 ™ 3 i n the second term and using equations (|25c|) , fl25d|) (only 



(ni+l)n 2 n 3 

a 



A (ni+l)n2ri3 
^0 



V^l + la" 



(A9) 



and 



(ni-l)n 2 n 3 



2 



(ni— I)ri2n3 



A 1 



nia . 



(A10) 



where 



^(ni+l)n 2 n 3 



fru; — (ni + 1) — n 2 hu 2 — n^huj^ — e 



and 



(ni-l)n 2 n 3 



ftcu — (ni — 1) 7tu;i — riihio^ — n^hus — ey. 



Replacing (|AD|) and ( |A10| ) in ( |A7| ) and doing the inverse Fourier transform for the am- 
plitudes with the same approximation as those made to get equation 
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G v = E i — (ni + 1) b?| 2 m \ P %\ 2 }. (All) 

1 l t a (ni+l)ri2n3 1 / l-fl I , . (m -I)ri2n3 1 lru J v ' 

n,occ 'i^O li/Ai/ 

If we do now a mean field approximation to take the phonon terms out of the sums, we 
get the following expression 

2 2 

G^ = (rn+l)V — r^-r, |Pv| 2 -nrE — IPof- (A12) 

n,occ 'iZAg n,occ llL\y 

The amplitudes and pj 00 are the only ones that contribute significantly to the first 
and second sum, respectively. Taking this into account and the fact that close to the 
resonance (ey — e ~ hui) A v ~ Aq 00 = A, equation ( |A12| ) can be written as 



of, = (w + 1) £ JL P y - wY, h: pI 00 ' ( A13 ) 



hA 



Finally the expression (|A13|) can be written as 



Gf, ~ w [ny-(ni + 1) -n^n\] = R w (A14) 
where iw = — - , 

occ occ 

The other terms appearing in ( [I5D can be transformed by the same procedure described 
above. 

In summary, to obtain the kinetic equations similar to that used in Q, i.e., involving only 
populations, we did not take into account the crossed terms. For example, if we consider 
the term with vj in (|25c| ) to obtain flAtp, then the expression for the amplitude a\ ni 1 



Jq 111 \[^Ql^\ ) IjU UUtdlll , lllCll U11C CApiCOSlUll 1U1 U11C ClllipilljUUC LlQ 

would have a term depending on vj . That term would describe a virtual process. In this 
process, an electron tunnel through the left barrier (j = 1) directly to the ground state 
(j = 0) emitting one LO phonon. Even when this process could occur, its probability is 
much smaller than that of the indirect process: first tunnel to the excited level, and then 
decay to the ground state emitting a LO phonon. Therefore, we make the sum of each 
process as being independent. The results obtained from equations flS^p are quite similar to 
those obtained from (|45|) . 
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FIGURES 

FIG. 1. Potential profile and energy levels at the resonant condition. bi(b r ) and d are the 
barriers and well width respectively. 

FIG. 2. The profile for the scattering region. 

FIG. 3. Diagram representing the electronic Hamiltonian (11). 

FIG. 4. Diagram representing equations (29). The vertical transition represents the LO phonon 
emission (absorption) . As the channels + and — have the same energy, they are represented in the 
same plane. The gray points represent the sites with the renormalized energies and the gray lines 
are the renormalized hoppings. 

FIG. 5. Diagram for the several rates appearing in (41a) and (41b). 

FIG. 6. Dashed line: the two polaronic peaks when 7 = 0. Solid line: if we assume 7 / 0, 
it appears a third peak corresponding mainly to one electron in the ground state with a pair 
(LO + TA) of phonons. 

FIG. 7. Solid line: the current for 7 / 0. For comparison, the electronic current when 7 = 
is shown in the inset (dashed line). It appears a shoulder in G (7 7^ 0) corresponding to the third 
peak of figure 6. The cross-sectional area used was S = 0.5 • 10~ 3 mm 2 . 

FIG. 8. The TA phonon population and the saser intensity S = ^3^3. 
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